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Abstract 


We study a class of VT -symmetric semiclassical Schrodinger oper¬ 
ators, which are perturbations of a selfadjoint one. Here, we treat the 
case where the unperturbed operator has a double-well potential. In 
the simple well case, two of the authors have proved in [B] that, when 
the potential is analytic, the eigenvalues stay real for a perturbation of 
size 0(1). We show here, in the double-well case, that the eigenvalues 
stay real only for exponentially small perturbations, then bifurcate into 
the complex domain when the perturbation increases and we get pre¬ 
cise asymptotic expansions. The proof uses complex WKB-analysis, 
leading to a fairly explicit quantization condition. 
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1 Introduction and results 


Operators that are 'PT-symmetric have been proposed in quantum mechan¬ 
ics as an alternative to selfadjoint ones. From a physicist point of view, it is 
of course very important to verify that the spectrum of such operators is real 
and there has been a considerable activity in this area [n, m, 0,0,0, 

cni, 0 , m, IS], m, ibi, m, m, m. m, m- . There is in particular 

an issue 0 of the Journal of Physics A which is devoted to non-selfadjoint 
operators in quantum physics, where the PT-symmetry property plays the 
main role. In a recent paper [Hj, E. Caliceti and S. Graffi study qualita¬ 
tively if, in a perturbative setting, the phenomenon called PP-symmetry 
phase transition occurs for PP-symmetric polynomial potentials, that is if 
the eigenvalues bifurcate from real to complex values. The reader may also 
find interesting references in that paper. 

It has been proved recently by two of the authors in 0 in the analytic 
category, that the PP-symmetric perturbation of a semiclassical Schrodinger 
operator with a real-valued single well potential, have real spectrum even 
for perturbations of size 0(1). Here we address the double-well case, where 
numerical results suggest that the situation is very different. For example, 
in Figure we have plotted the eigenvalues of the Schrodinger operator 
Pe,h = h‘^D^ + 0.05x^ — .5x^ + iex in a neighborhood of the barrier top, for 
h = .01 and e = A:.10“”*, k = 1,2, 3,4 and m = 2,3,4,5. They have been 
computed by a simple finite difference scheme using the package scipy in 
python, and plotted using the package matplotlib [Hj. It appears that the 
eigenvalues close to Pq > 0 are real, giving a numerical illustration of the 
results in 0 since we are then in a simple well situation. In the double well 
case, that is for Eq < 0, the eigenvalues close to Eq seem to be non-real 
except for extremely small e, and this is the phenomenon we consider in this 
paper. 

We recall that an operator is said to be PP-symmetric when it commutes 
with the operator PP, where P is the parity operator, and P the time- 
reversal operator given by 

(1.1) Vu{x) = u(—x) and Tu(x) = u{x). 

Here, we study small perturbations = Pe(x,hD) of self-adjoint semi- 
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^ xin-^ £ = k.lO TO = 6 



- 0.1 0.0 0.1 
Re(A) 


^ xlQ-^ g = ^'-lO m = 3 



-1 -'- 

- 0.1 0.0 0.1 

Re(A) 


^ xlO“® ^ m = 9 



-1 -^- 

- 0.1 0.0 0.1 

Re(A) 


Figure 1: Eigenvalues of + O.OSx"^ — .5x^ + iex near 0, with 

h = 10“^ and e = A:.10“"*, A = 1,..., 5. 


classical Schrodinger operators of the form 

(]2 

(1.2) P,{x,hD) = -h^^ + Ve{x), Ve{x) = Vo{x) + ieW{x), 

where Vq and W are smooth functions on M and Eg G M a fixed energy, 
satisfying the following assumptions: 

(Al) Vq is C°° and real-valued on M. 

(A2) There exists uiq > 0 such that, with {x) = (1 -|- x^) 2 , 

VA E n,3Ck > 0, Vx E M, |eJ^^(x)| < Cfc(x)”*°“^ 
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and 


VxGM\]-Co,Co[, \Vo{x)-Eo\ > 

(A3) For some Eq G M, the equation Vo{x) = Eq has exactly four solutions 
ai < l3e < Pr < ar, with V^iae) < 0, Kil^r) < 0, and 

Vo'(ar) > 0. 


Therefore, the operator Ph^ on L^(M) with domain 
V = {uG {x)^°u G t2(M)}, 

is a self-adjoint Schrodinger operator with a double-well potential, and its 
spectrum consists only in real eigenvalues with multiplicity 1. 



Vo{x) 

Eo 


Figure 2: The unperturbed potential Vq 


Concerning the perturbation isW, we suppose that, with mo given in 
(A2), 

(A4) VF is a C°°, real-valued function on M, such that 

yk G n,3Ck > 0 , Vx G M, < Ck{x)^°-^. 


Then, for any e g] — eO) £o[ with eo > 0 small enough, the spectrum of the 
unbounded, closed operator Ph^e on L^(M) with domain V, is still discrete 
in a complex (h-independent) neighborhood of Eq. We are interested in the 
semiclassical asymptotics (i.e. as h ^ 0) of the eigenvalues of Ph,e near Eq, 
and we want to know whether they stay real or not for e > 0 small, under 
the following PT-symmetry assumption: 
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(A5) The operator is PT-symmetric, that is 


[Ph,e,rT]=0, 


where V is the parity operator, and T the time-reversal operator given 


in (1.1). 


In terms of the real and imaginary part of the potential V^, this assump¬ 
tion is equivalent to the property that Vq is even and W is odd. In particu¬ 
lar, when e = 0 we have a symmetric double-well potential, and the turning 
points satisfy 


(1.3) 


ai = -ar, /3i = -/3r. 


Since we are going to use complex WKB constructions in a neighborhood 
of the wells, we need to suppose also that 

(A6) Vq and W have analytic extensions to a neighboorhood [/ in C of the 
convex hull of {x G M; Vo{x) < Eq}. 

Remark 1.1 In (A2), and the corresponding assumption on W, we can also 
treat the case mo = 0. Then Phfl is still self-adjoint, and its spectrum in 
] — oo,Eq + 1/C[ consists also only in eigenvalues. Our results remain valid 
in this case. 

In this paper, we will first obtain a Bohr-Sommerfeld like quantization 
condition for the eigenvalues of Ph,e in E{Eo,l/C). This result does not 
rely on the "PT-symmetry assumption on Ph,e- 

Let us define the relevant action integrals. By the implicit function the¬ 
orem, there exists £o > 0 such that the equation V^(z) — E = 0 has exactly 
four solutions for {E,e) G D{Eo,eo) x Zl(0,eo)) that we denote ai{E,s), 
(3i{E,s), j3r{E,e) and ar{E,e), depending analytically on E,£, in such a 
way that. 


(1.4) 


a,{Eo,0) = a,, l3,{Eo,0) = P,, • = £,r. 
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These analytic functions of {E, e) are called turning points at energy E. The 
action integrals are the functions given by 

rME,e) J 

(1.5) I.{E,e)= {E-Vs{t))- 2 dt, 

J a,(E,£) 

where • = £,r and the branch of 1 1 —)■ (£' — 14 (f ))2 is the one that is real and 
positive for {E,£) = {Eq,0) on the segment between a, and /?,, • = £,r, 
and 

rdr{E,£) ^ 

(1.6) J{E,e)= {Vsit)-E)2dt, 

Jdi{E,e) 

where t i—)■ (14(t) — E) 2 is real and positive for {E,£) = {Eq,0) on the 
segment between /?£ and jSr- These functions Ig, Ir and J are analytic with 
respect to E, s. 

We shall say that a function f{z,h) defined in nx]0,/io], analytic with 
respect to z G a domain of C*^, has an asymptotic expansion 

j>o 

in Hol(n), when the /j’s are holomorphic in and, for all compact subsets 
K cn and all iV G N, 

N 

(1.7) \\f{z,h) = 0{h^^^). 

j=0 

For a function k{E,e) we shall denote, for e G M, 

(1.8) k*{E,e) = k(E, -e) and k\E,£) = k(E,s). 


Theorem 1.2 Assume (A1)-(A4) and (A6). There exists a fixed neighbor¬ 
hood = Q 1 XQ 2 of {Eq, 0) in CxC and holomorphic functions 7^(1?, e; h), 
75-(FI, £', h) with complete asymptotic expansions in Hoi (fl) and of the form 
1 + 0{h) such that if 


(1.9) 


f{E,s;h) 


1 ^ 

— sin (/^//i) sin(/r//i), 


^-ilvlh^r 
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then when h > 0 is small enough and e G fii n M, the eigenvalues of 
in coincide with the zeros of f{-,e]h) in the same set. The algebraic 
multiplicity of each such eigenvalue E coincides with the multiplicity of E 
as a zero of f {■,£', h). 

We may assume that ri is invariant under the map {E,e) i—)• {E,e), and 
we have 

(1.10) 1} = k, i; = Ir, r = J, ( 7 I)* = li, {Y-Y = 


Assuming also (A5), we have 

(1.11) /{E, e) = pp^ cos{I/h) cos(lV^) — sin(///i) sin(I'^//i), 

where p{E, e; h), I{E, e; h) have complete asymptotic expansions in Hoi (H) 
with p = l + 0{h), I = I + 0{hY. Further, I* = I, p* = p and J* = J = J'^. 


Remark 1.3 The function f in (1.9) can also he written 


(1.12) f{E,e]h) = p(,prCos{Ii/h)cos{Ir/h) — -e sm{I(^/h)sm{Ir/h), 


where 

h — h, Ir — Ir = Pi — 1, Pr — 1 = 0{h). 

These quantities have complete asymptotic expansions in powers of h and 
enjoy the symmetries, 


p} = pe, p; = Pr, /; = le, i; = u. 


U nder the additional assumption (A5), we have pr = p}, E = 1} and we 
get (1.11) with p = pi, I = E, I = E. 


The zeros of E cos{E/h) and E 1 —)• cos{Ir/h) can be thought of as 
“approximate eigenvalues for the left and right potential wells respectively”. 
They are given by the Bohr-Sommerfeld quantization conditions, 

(1.13) 2E{E,e;h) = {2k + l)7rh, A: G Z, 
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and 


(1.14) 


2Ir{E, e; h) = {2k + l)irh, A: G Z, 


respectively. The last term in (1.12) is exponentially small and represents 
the tunnelii ig int er actio n between the potential wells. The solutions E = 
E*{£-, h) of (1.13), (1.14) are situated on the curves r,(e; h) dehned by 


(1.15) 


lmI,{E, e; h) = 0, 


respectively, and the distance between consecutive zeros E* and is of 
the order h. Here, we use that, for • = i,r, 


BeI, 


1 

2 



{E -Ve{x))-^dxi^Q. 


When e = 0, T, are real neighborhoods of Eq and the E* are real. We notice 
that 


ReEi{e) = Ei{0) + O{e^). 


It is clear that the set of zeros of f{-,e;h) is (in a suitable sense) expo¬ 
nentially close to the union of the solutions of ( 1.13| ) and (1.14). Below we 
give a such a detailed result in the (PT-symmetric case. 


Now we adopt the assumptions (A1)-(A6) as well as the following as¬ 
sumption: 


(A7) We have 


Notice here that 



{Eo -Vo{x))-^W{x)dx ^ 


0 . 


fdi 1 

2ideIe{E,£)= {E-Ve{x))-2W{x)dx. 

J on 

Possibly after the substitution (e, IP) i—)■ (—e, —IP), which does not change 
Pg, we may assume that the integral in (A7) is > 0. 

Under the assumption (A7), we have = Ir and hence for real e, that 
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Theorem 1.4 We make the assumptions (Al) to (A7). The values Ek ;= 

Ef. and E^ = Ek are situated on the curves F = and F = F^, where F is 
of the form 

(1.16) ImE = g(ReE,e]h). 

Here 

(1.17) 9 {t,e; h) ~ g{t,£) + hgi{t,£) + ... in Hoi (neigh (E^o, 0), C^) 
is real for {t, s) real, and we have, 

(1.18) g{t, £) = 0) + 0{h‘^)^ £ + O(e^). 

There exists a fixed neighborhood H = Hi x H 2 of {Eq,0) G C x M, such 
that for £ £ 0,2 and for h > 0 small enough, 

(1.19) f-\0,£) c[Jd (^Ek,r{Ek,£)) u\Jd (Wk,r (Wk,£)) , 

k k 

where r{E,£) < (7/je“ is given by 

(1.20) r{E,£) = Chmin {l,max{h/£, T)e-'^AE,e)/hi^ ^-Kej{E,e)/h^ 
and C > 0 is large enough. Moreover, 

— when these discs are disjoint, f{-,£) has precisely one zero in each of 
D (^Ek,r (^Ek,£^^ and D (^Ek,r (^Ek,£^y 

— in general /(•,£) has precisely 2 zeros in 

D (Ek,r (^Ek,£^^ U D (^Ek,r (^Ek,£^^ . 


We finally discuss the more precise behaviour of the eigenvalues when |e 
is exponentially small. The function 


f/{pp^) = cos 




1 


2J/h 


sm 
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is real-valued on the real axis and has a sequence of local minima Ek{£', h) 
such that 


^fc(0; h) = Ei{Q; h) + O ^ 


where E^(Q] h) = E'^{0; h) is defined in (1.13), (1.14) and we know that when 
|e| < there are two eigenvalues of Ps exponentially close to i?fc(0; h) 

and that we obtain in this way all the eigenvalues in a fixed neighborhood 
of Eq. It will be convenient to introduce a “Floquet parameter” k G M and 
to set 


f{E,£,K;h) = 




We still have a sequence of local minima Ek{£,K;h) satisfying Ek{£,K -|- 
Tr;h) = Ek_^_l{e, k; h) and the zeros of f{-,£,K;h) are now confined, two by 
two, to exponentially small neighborhoods of Ek{£, n] h). We concentrate 
on one such local minimum Ec{£,K‘,h) = E}f{£,K]h) and we restrict the 
attention to a “window” 


(1.21) E = Ei + hF, £ = h£, K = K + I{Ei,0)/h, 

where Ei is a real parameter G neigh (Flo j 1^)- Assume that Ec{£,K-,h) be¬ 
longs to the window, so that 

(1.22) Ec{£, K] h) = El + hFc{£,K,; h), 

where Fc is the corresponding critical point of / in the variables F with 
£, K as the new parameters. Thanks to the rescaling, this critical point is 
uniformly nondegenerate. / and Fc are even functions of e and it follows 
that 

(1.23) Fc{£,K-h) = Fc{0,F,h) + Oi^). 


We make the assumptions (Al)-(A7) and discuss the zeros of / near 
the critical value Ec{£,K]h) = Ei + hFc{£,'ii\ h). The key point is that, in 
this regime, we are able to write the quantization condition as a second 
order polynomial (up to a non-vanishing factor), with a sharp control on 
the coefficients. 
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Theorem 1.5 The critical value f^{£, k; h) = f{Fc, e, k; h) is of the form 

(1.24) = m{£,K;h){^ -£c(K]hf), 
where 

(1.25) ec = ^(k; 

Here, £, m are classical symbols of order 0 as in 
satisfying 

(1.26) ^ 2j57(k~0)I’ 

Further, 

(1.27) f(F,£,K;h} = f^(£,K;h} q(F,£,K; h)(F - Fc(£,k, 1; h}f, 

where q is a symbol of order 0 and 

(1.28) q{Fc,0,li-,0) = 2idEl{E,i0,K,0),0)f. 

/(.,£, K;/i) has two zeros in a small neighborhood of Fc when counted 
with their multiplicity, and 


1.7), with leading terms 


— when |el < £c{ii; h) the zeros are real and simple, given by 

(1.29) q{F,£,K;h)^{F - Fc{£,K,l]h)) = ±{-^{£,71; h))^ 

— when |el = £c(k; h) we have a double zero, 

(1.30) F = Fc. 

— when |e| > £c{k‘, h) the zeros are simple, non-real and complex conju¬ 
gate to each other, given by 

(1.31) q{F,£,K-h)^{F - Fc{£,K,T,h)) = ±i{f{£,K-,h))^. 


When K = K + I{Ei,0)/h belongs to vrZ, these values give the eigenvalues 
of Pe near Ec{£, k; h) via {1.21). 


Eventually, we would like to mention the paper m by C. Gerard and A. 
Grigis, where the authors study the eigenvalues of self-adjoint Schrodinger 
operators with a double well potential. They also obtain a quantization 
condition, using what they call the ’’exact WKB method”. Our method 
here is slightly different and more explicit about the connection formulas at 
the turning points. 
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2 The complex WKB method 


We recall here briefly elements of the complex WKB method in a general 
setting. Consider a Schrodinger equation, 

(2.1) —h?u"{z,h) + V{z)u{z,h) = ^, 

in a bounded, simply connected open set 17 C C where the potential V is 
holomorphic. We look for a solution of the type 

(2.2) u(z,h) = a(z,/i)e*^(")/'^, 

where a{z, h) has a formal asymptotic expansion in a sense to be defined 
later on, 

H-OO 

(2.3) a{z, h) ^ aj{z)h^, 

j=0 

and the Oj’s are holomorphic functions in U. The function (p is called the 
phase of the solution u, and a{z, h) is called its symbol. 


A function u{z, h) of the form (2.2) is a solution to (2.1) if and only if 

(2.4) + K(z))(e*^W/^a(z, h)) = 0, 
or 

(2.5) {—{hdz)‘^ — 2ip'{z)hdz — ihp"{z) + p {z)"^ + V{z))a{z, h) = 0. 

If is a solution of the eikonal equation 

(2.6) ip'{zf + V{z) = {), 


then (2.4) is equivalent to 

(2.7) 


'f'{z)dz + - y a{z, h) = 0. 


Replacing a{z,h) by its formal asymptotic expansion (2.3), and canceling 
successively the powers of /i, we obtain a sequence of transport equations 


( 2 . 8 ) 


(p'{z)dz + ^(p''{z)j ao = 0 , 

^\z)dz + \ 3 > 1 - 
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Definition 2.1 A formal WKB solution u^kb of the equation (2.1) in U is 


a pair {(p, (oj)) of an analytic function (p in U verifying the eikonal equation 
(2.6), and of a sequence (aj) of analytic functions in U which satisfies the 
transport equations 1(2.8\). We denote it 


(2.9) 




j>0 


We suppose from now on that V{z) 7 ^ 0 for all z € U. Then we fix a 
determination of z (z))^ in U, and we solve the eikonal and transport 

equations in U. 


Proposition 2.2 The solutions of the eikonal equation (2.6) are analytic 
functions in U, and they can be written 


( 2 . 10 ) 

for some zq G U, and some C G C. 


(p{z) = ± / {—V{w))^dw + C 
J ZO 


Now we fix a such a solution ip in U. It is then easy to prove by induction 
that, given an initial data, the system of transport equations has a unique 
solution. Therefore we have the 


Proposition 2.3 Let (a^ )^Q he any sequence of complex numbers. Then 
the Schrodinger equation (2.1) has a unique formal WKB solution in U, such 
that 

Mj G N, aj{zo) = 0°. 


Moreover, the function oq is given in U, for some suitable constant C G C, 
by 

ao{z) = C{ip\z))~^/‘^. 


We want now to associate true solutions of the Schrodinger equation (2.1) 
to the formal ones we have constructed above. It is convenient to introduce 
the notion of Stokes line for the potential V. 


14 







Definition 2.4 Let U he a simply connected open set in C where V is 
holomorphic. A Stokes line is a curve a : I ^ U such that 

/ t ^ 

{—V {a{T)))'^ a' {T)dT = 0, 

for all s,t G I- Here, I is any interval starting at 0 and ending at 1. 


Notice that 

where is a solution of the eikonal equation. Thus, a Stokes line is nothing 
else than a level curve in U of the imaginary part of the phase g). 

The following proposition is well known (see for example [23] for a proof), 
and can be considered as the fundamental rule of the complex WKB method: 
always move in a direction where the modulus of the phase factor increases, 
thus in particular transversely to the Stokes lines. 


Proposition 2.5 Let U he a simply connected bounded open subset of C, 
such that V{z) ^ h for all z G U. Let ip be a solution of the eikonal equation 
in U. Let also j :]0,1[—)■ [/ be a curve in U such that 

(2.11) Vt g]0, 1[, ^(-Imv7(7(t))) > 0, 


Then there exists a neighborhood LI G U of y such that, for any formal 
WKB solution Uwkb{z,h) = ’ there exists a solution u 


of the Schrodinger equation [2.1) in such that 


( 2 . 12 ) 


u(z,/i) = 


where a is holomorphic with respect to z G Ll, and 

(2.13) a{z,h) ~ '^^aj{z)h^ mHol(n). 

k>0 
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3 WKB analysis near a simple turning point 


In this section we follow closely the presentation in [23] . Let II C C be open 
and simply connected, and V £ Hol(n). We suppose that V has a unique 
zero zq in O, and that it is a simple one: 

(3.1) V{zo) = 0, V\zo) ^ 0. 


We are interested in solutions u in 0 of the general Schrodinger equation 

u{z,h) = a{z,h)e^^^^/’^. 


(2.1) of the form 


(3.2) 


Notice that, contrary to (2.2), here we have chosen not to put the factor i 
in the exponent to simplify the notations. For the same reason, we will also 
assume that zq = 0. 


As in Section]^ we obtain first the eikonal equation 
(3.3) ^Xz) = V{z)K 

in H. By assumption there exists a function F, holomorphic in 11, such that 

V(z) = zF{z) 

and F{z) ^ 0 for z £ Vl (we may decrease II whenever necessary). It is 
therefore clear that ^{z) is multi-valued in general, and to better understand 
the structure of this singularity we pass to the double covering H of H \ {0}, 
setting z = w'^. Then 

d Id 
dz 2w dw ’ 

and if we set _ 

r V{w) = V{z) = ViP‘F{vLp‘)^ 

i ^{w) = (p{z), 
the eikonal equation becomes 

dw^{w) = 2w‘^F{w^)'2 . 

Notice that the right hand side is an even holomorphic function. If we also 
require that <^(0) = i^(0) = 0, we see that is an odd holomorphic 

function of the form 

2 11 
<^{w) = -F{nP‘)w^, where F{0) = F( 0)2 = I/'( 0 ) 2 . 

3 
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In the original coordinates, we get the double-valued solution, 


(3.4) 


2 ~ 3 

<^(^) = ^F{z)z2. 


Now we study the Stokes and anti-Stokes lines having 0 as a limit point. 
Since, with respect to Sectionj^ we have removed the factor i in the exponent 
in (3.2), Stokes lines are now level curves of the real part of ip, and level 


curves of Im p are called anti-Stokes lines. On such curves we have Re (/J = 0 
or Imp = 0, which is equivalent to \mp^ = 0, and to \mz‘^F{z)‘^ = 0. In 
other words, these curves are given by 


(3.5) 


{z G n, 3t £ R, z^F{zf‘ = t^}. 


Taking the cubic root, we see that Stokes and anti-Stokes lines reaching 0 
in the limit, are contained in three curves 7 ^ given by 

(3.6) 7 fc = {z G 12 , 3t G M, zF{z)^ = k G {0,1, 2} ~ 'Ll?,!.. 

In the case where 14^(0) > 0, the situation is as shown in Figure]^ Each 
curve 7 fc \ {0} is divided into a Stokes line 7 ^ (plain lines) and an anti-Stokes 
line 7 ^ (dashed lines). The three Stokes lines delimit three closed Stokes 
sectors k G Z/3Z, where is the sector that contains 7 ^. In Figure]^ 
we have also drawn a Stokes line inside each sector. 



Figure 3: Stokes lines close to a simple turning point 


For k G Z/3Z, we denote by pk the branch of <y 9 in 12 \ 7 ^ such that 
Pk{^) = 0, and Re^ifc < 0 in Notice that pk+i and pk are both well 
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defined in U Sfc+i and satisfy 


(3.7) 


ipk-\-\ — in S/j U 


According to Proposition 


Schrodinger equation (2.1) in such that, in 


2.5, there are solutions u 


Uk, k G Z/3Z, of the 


Uk{z,h) = akiz,h)e'^^^''^/’^, 

ak{z,h) ~ '^akj{z)h^ in Hol(i;fc). 
j>o 


This asymptotic description extends to 11^, the complement of an arbitrarily 
small neighborhood of 7 ^U{ 0 }, that can be reached from by crossing the 
Stokes lines transversally. We also recall that akfi is unique up to a constant 
factor and that we can choose 

(3.9) ak,o{z) = Mz))-'^, 


for any branch of the square root. 


Recall that if u, v are solutions to the Schrodinger equation, then the 
Wronskian 

Wh{u,v) = {hdzu)v — u{hdzv). 

is constant, and vanishes precisely when u, v are collinear. Let j,k G Z/3Z. 
Applying the asymptotics of Uj and Uk at some point in the interior of 
SjUSfc, we see that, recalling (3.9), Wh{uj,Uk) has an asymptotic expansion 
in powers of h, whose first term is given by 


(3.10) Wh{uj,Uk) = 2ajflakflip'j + 0{h) = 








+ 0{h). 


We fix a branch of {^'k)'^ in '^k for each k G Z/3Z. For two different Stokes 
sectors, Sj 7 ^ we have in the interior of Sj U that 


(3.11) 




for some Vj^k £ Z/4Z which are odd, and such that Vj^k = ~^k,j- Thus, 
starting from Sq, we can make a tour around 0 in the positive direction and 
we get that 

in El, 

(3.12) ((^' 2)2 = in S 2 , 

(V^o)^ = (</^2)^ in So. 
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This means that if we follow a continuous branch of 2 around 0 in the 
positive direction, then after a turn, we obtain the new branch 


(3.13) 


(l^O, 2 + 1^2, l+! 2 l,o) 




But (</?o )2 = for a suitable branch of the fourth root, and if one follows 
this function around 0 once in the positive direction, we obtain This 

gives the co-cycle condition 


(3.14) 


vqp + 1 ^ 2,1 + 1 ^ 1,0 = -1 mod 4. 


(3.10) and (3.11) imply that 

Wh{uj,Uk) = 2i'^^’>^ +0{h). 


Now we describe the linear space of solutions of the Schrodinger equation 


(2.1) in n. It is of course of dimension 2, and any two of tt_i, uq, ui are 
linearly independent, so we have a relation 


(3.15) a_iu_i + aouo + aiui = 0, 

where the vector (a_i,a 0 )«i)^ G \ {0} is well defined up to a scalar 
factor. Applying W{uj, •) to this relation, we get 


(3.16) 

(^tlj , Uk)) j^k j 

'a-i\ 

ao 

which is 

\ 

a system of the form 

^ «i / 


/ 0 a b\ 

/a_i^ 

(3.17) 

1 —a 0 c I 

ao 


\-b -c ol 

V ai J 


= 0 , 


= 0 . 


The triplet (a-i, ao, ai) = (c, —b, a) is a solution, so up to a common factor, 
we have 


(3.18) 


a. 


= ±i + 0{h). 


More precisely, the values of a, b and c are given by the equation (3.15), and 
we get, after inserting a factor 1/2, 


(3.19) 



/ \ 

+o{h) 

\ / 


/ \ 

+ 0{h). 
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Rema rk 3 .1 Sometimes it is more natural to change the notation, writing 
i(pj in (3.8) instead of (pj so that Uj{z] h) = aj{z-, with Im ( p,- > 0 

in Yij. Then (3.9) becomes ajfl{z) = = V{z)~^/^ and in (3.11), 

(3.12), if'- must be replaced by itp',. 


4 WKB solutions near the wells 


From now on, we consider the equation Ph,eU = Eu, that is 
(4.1) -h^u” + {Ve{x) - E)u = 0, 

where 14 = Vq + ieW satisfies (Al) to (A4) and (A6). For the moment we 
do not assume the PT-symmetry property (A5). 


Let us now define some formal WKB solutions to the Schrodinger equation 


(4.1) near the wells. For {E,e) G D{EQ,eo) x 11(0, eo)) the equation 14(x) = 
E has exactly four solutions in U, the domain of holomorphy of 14, that are 
called turning points at energy E. We have denoted them ai{E, s), ^e{E, e), 
j3r{E,e) and ar{E,e), with. 


(4.2) a,{Eo,0) = a,, ^.(Kq, 0) = /?,,• = £, r. 

We have drawn in Figurej^a typical configuration of the Stokes lines starting 
at each of the turning points, when E ^ Eq and £ ^ 0. 


We shall work in the cut complex plane along \ai,j3(\ U or more 

precisely in the cut version U of U, so that we have two determinations of 
X I-)- (14 (x) — E) 2 in IJ. We denote 

X (14(x) - K)!, (resp. x ^ (14(x) - £')4, x ^ (14(x) - E)?), 

the determination which is real and positive for e = 0, E = Eq and x G 
] — (X),a(\ (resp. x G]/?^,/?^), x G]ar,+c>o[). Notice that 

Vx G tj, (K(x) - E)j = (K(x) - E)} = -{V,{x) - E)i. 


First, we concentrate on the situation near the left well. We denote Eg, 
Ef and the three Stokes sectors near ai, and Sq, S'f and those near 
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Figure 4: Stokes lines and Stokes sectors 


/?£. For eac h sector (resp. 5^), k G Z/3Z, we choose a solution (resp. 
of (4.1) such that 


(4.3) 


ui{z, E, e, h) = ai{z, E, e, h)Eviiz,E,e)/h 
vi{z,E,e,h) = in 5^ 


Here (resp. V'D is a solution of the eikonal equation 
(4.4) = V,{x)-E, 


vanishing at z 
such that 

(4.5) 


ai{E, e) (resp. at /3i{E, e)) for {E, e) G D{Eq, e) x 11(0, eo)) 

Vz G Re{iip[{z,E,£)) < 0, 
yzeSl Ke{i^i{z,E,e)) <^. 


The amplitudes a|, and b\, in (4.3) have asymptotic expansions in Hol(S^) 
and Hoi (5^) respectively, in the sense of (1.7) The phase functions ^p\. and 
are well defined in and V’^+i are well defined in 

where they satisfy 


(4-6) (pi — —ipi — —V’fc+i- 

We choose the functions (resp. f^) so that they are holomorphic also 
with respect to {E,£) in D{Eq,£q) x 11(0, eo). The asymptotic expansions 
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of the amplitudes and 6 ^ extend to x D{Eq,£q) x L)(0,eo) and x 
D{Eq,£q) X L)(0,eo) respectively, where (resp. is the complement of 
an arbitrarily small neighborhood of 7 “^ ^ (resp. 7 ^^ in [/ (see Figure j^. 

We fix now a choice for the principal symbols q and h^. g of these six 
solutions and k G Z/3Z. 


By (|4.4|) and (|4.5| we have 

i(pQ{x,E,£)= f {Ve{t) - E)ldt, 


(4.7) 




^{x,E,e) = - r{Ve{t)-E)jdt, 

J OLp 
pX 

iip[{x,E,£) = - (Veit) - E)ldt. 

. J at 

Then we choose the principal symbol Og g of Og to be 


(4.8) 


a^o,o(x,i?,e) = [(i/o)']"" = (K(x) - i?), F 


In order to fix the principal symbol of and we have to choose 

1 ^ 0 , 1 ) 2 ^ 1,-1 and 2 ^- 1,0 in Z/4Z, odd, such that (3.14) holds. We take 

(4.9) z/ 0,1 = 1, 2 ^ 1-1 = -1, and = 1- 


Then by (3.12) we have. 


1 


(4.10) 0-1,o(a;, e) = e) = -*(i4(x) - .F)^ 

and, 

(4.11) af o(x, E, £) = g(x, E, e) = i(K(x) -E)^V 

With these choices, we have 

(4.12) Uq = r+(/i)til;^ + T_ih)u{, 


for some symbols T±{h) with T±{h) = 1 + 0{h). Without changing the 
leading asymptotics, we can replace by and we get 

(4.13) Uo = '*^-i+^i- 
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Let us now consider the solutions v^, k € near /?£. By (4.4) and 

(4.5) we have 


(4.14) 


i^pQ{x,E,e) = 

i'il;{{x,E,£) = 


- / (Veit) - E)ldt, 

Jfh 

-- [ (Veit) - E)idt, 

Jdi 

[ {Ve{t) - E)idt. 


We here consider — z as the basic variable for the Schrodinger equation and 
correspondingly, we choose the principal symbol 5 q g of 6 q to be 

(4.15) bloix,E,e) = [(-i^^)']-5 = {Ve{x) - E)-\ 

and we fix the principal symbol of vf and choosing and 

as in ( |4.9[ ): 

(4.16) z/ 0,1 = 1, = -1, and = 1- 

We get, 

(4.17) b\oix,E,s) = ^6^g g(x,^,e) = -i(W(x) - 77)“^ 
and 

(4.18) bl,{x,E,e) = E,e) = iiVs{x) - E)-\ 


We further fix a choice of The principle of the WKB method ensures 
that we can choose v{ to be proportional to and v^_i to be proportional 
to u\. Notice first that 

i'il>{{x,E,e) - iiptiix,E,£) = ( {V^{t) - E)^dt =-iIi{E,h), 

J a.^ 

where we have set 

(4.19) h{E,e) = f\E-Ve{t))idt 

Joi( 
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where {E — V£{t))w is real and positive for e = 0, i? G M close to Eq and 
t G]a£,/3f[. The same way, we see that 

ii’-iix, E, e) — i(p\{x, E, e) = iIi{E, h). 


Concerning the principal symbols, we have first 


(4.20) 


{V£{x)-E)j^ = -i{V£{x)-E) 

{V£{x) - E)j_ = iiVeix) - E)i 


1 

4 

m ? 


_i _i 

Here we have denoted {Ve{x)—E)j^_^ (resp. (V):(x)—) the determination 

1 — — ~ 

of {Ve{x) —E)~i obtained on ]/3i, /3r[ by extending {Ve{x) — E)^ on [/ along 

a path in the upper (resp. lower) half plane. Thus we see that 
(4.21) 6^ g = and feli^g = —g. 


and we can assume that 




u{ = 


for some symbols cr^ such that = 1 + 0{h). After replacing by 
(which does not change the leading asymptotics) we may assume 
that cr^ = 1: 


(4.22) 


u\ = 


The same discussion applies to the solutions associated to the well to the 
right. The main rule is simply that the right well becomes the left well of 
the operator = —h?d‘~ + 14(—5?) under the change of variables x = —x, 
and we let u^., be obtained from the corresponding null solutions vl 
of Pg — E. Note that the Stokes sectors SI correspond to the sectors 
S|, si to the left, defined exactly as Sl- (Cf. Figure |^) 


This means that we have the 6 solutions u^, k which satisfy (4.3) 
with replaced by “r” and V’fe solutions to the eikonal equation 
(4.4), vanishing at z = ar{E,e), z = f5r{E,e) respectively, satisfying (4.5) 
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4.8), 


with ‘T’ replaced by “r”. The principal parts o^q, 6^ q are given by 

dilol), prTl), pT5l), (|4T7l), dllsl) with “f’ replaced by “r”. Here {V^ x)- 

— ^ - - - - 

E)r is the branch which is > 0 to the right of when e = 0 and E is real 
and close to Eq. Again, we can modify the choice of by constant factors 
1 + 0{h), so that the analogue of (4.13) holds: 

(4.23) u^o = uLi+u[. 


Then we can modify by constant factors 1 + 0{h) so that (4.22) holds 
with ‘T” replaced by “r”: 


(4.24) 


uLi = -iEE{E,s)/h^r^ 
Ul = ie-iIr{E,e)/h^r_^^ 


where now (cf. (4.19) 
(4.25) 


IriE,e) = 


par 

/ {E 

J 8r 


- Ve{t))i,dt 


and {E — V^)w is now defined near the right well ]j3r,oir[ as the branch of 
the square root which is positive on ]/5r,ar[ when e = 0 and E is real and 
close to Eq. 

Notice that this fits with the principle of transforming everything from 
the right to the left by putting 14 (x) = 14 (—x). 14 has the left turning 
points ai = —ar, Pi = —/3r and 

[^\e - V,ix))Ux = - f^\E - 14 (x))idx = Ir, 

J aj^ J ar 

with the natural branches of the square root. 

Now we build two convenient independent formal WKB solutions Wq and 
Wq near the barrier. We set 


(4.26) 

Since 




vf = i(V, - + 0(h)) = i(V, - E)m^e-^^o/h(l + 0(h)), 
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and 


/_ 1 = -i{V, - + 0{h)) = -i{V, - E)m^e-^^o/\i + 0{h)), 

we have 

On the other hand, we have 


r 

^0 


= {V,-E)m^E^y\l + 0{h)), 


wi = 


so that 

where Si^r = 1 + 0 (/i) is a symbol, and 

(4.27) J{E,e) = -ii;l-ii;i= / {Ve{t) - Ey^dt. 

jpt 

We now replace Vq by di^rVg (which does not modify the leading asymptotics) 
so that 


(4.28) 




Notice that (4.28) fixes a choice for the formal WKB solution Vq. 

In the same way, we set 

(4.29) 

and we have 


< = 4(t.: - oil). 


Wo = 

then replace Vq by dr/v^ and get 

(4.30) 

which we take as the definition of Vq. 


In analogy with the equation prior to (4.11) we have 

(4.31) v^o = lUh)vi+7iih)vy, 


(4.32) v^, = yyh)vy+Y_{h)vi 

where 7 ^ = l + 0{h). Having already adjusted by factors l + 0{h), there 
is no place for further adjustments, so we have to refrain from the possibility 
of replacing 7 ^ by 1 . 
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5 WKB expansions of L?‘ solutions outside the well 


In this section we focus on solutions of the Schrodinger equation (4.1) that 
are in a neighborhood of +oo or —oo respectively. The existence of such 
solutions follows from the general theory of partial differential equations 
since Ph^e — E is elliptic at infinity for E G D{Eq,£q). We are interested in 
their asymptotic behavior as /i —)• 0. 


Let 6o > 0. We consider the eikonal (2.6) and transport equations (2.7) 


on the half-line \ar -|- (5o, -|-oo[. There exists eo = eo(5o) small enough, such 
that for all {E^e) G D{Eq,£q) x D{0,£o), and for x G]a,. -|-(5o,+oo[, the 
function 


(5.1) 


ip+ix,E) = i 


J a 


ar{E,£) 


{Ve{t)-E)}dt, 


is a smooth function of x, and an analytic function of {E,£), which solves 

1 

the eikonal equation. We recall that t i—>• {Vs{t) — E)r is real and positive 
in ]ar + 6o, -|-oo[ when e = 0, so that 


(5.2) 


Re(i(^+(x, E)) < 0. 


It is then straightforward to obtain the existence of the solutions of the 
corresponding transport equations, and we get the 


Proposition 5.1 For all do > 0, there exists eo > 0 such that, for all 
{E, e) G D{Eq, £q) X E(0, eo), the equation (4.1) has a formal WKB solution 
u+ in ]ar + do, +oo[, 


(5.3) 


u+{x,£,E,h) = e*‘^+(^’®)/'^^aj(x,L;,e)/i^ 

j>o 


where aj is with respect to x G]a 2 + do, +oo[, and is holomorphic with 
respect to {E, e) G D{Eq, eo) x 11(0, eo). Moreover we can choose uq so that, 


(5.4) 


ao{x,E,e) = 


and we have the estimates, for all j > 0 and all k > 0, 
(5.5) |af^(x,L;,e)| = C>((x)“^“^'^^+^"*-^)“^). 
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Proof: We prove (5.5). Differentiating the eikonal equation 


ip'+{xf = E- Ve{x), 

and using (A2), we easily obtain by induction that, for A; > 1, 

(5.6) = 

Of course we also have, by (A4), 


so that (5.4) gives 

(5.7) 


(x) ° <\ao{x,E,£)\<{x) 4° 


Then, differentiating the hrst transport equation 

(5.8) 


m" fx) 

:p'+ix)aQ{x,E,£) + + aoix,E,£) = 0 , 


we obtain by induction that, uniformly for {E,£) G D{Eq,£q) x 74(0, eo), 

(5.9) |af (x,ii;,e)| 


which proves the estimates for j = 0 and all A; G N. Now suppose that (5.5) 
holds for some j and all A: G N. The {j + l)-st transport equation is. 


and, setting 

(5.10) 




Oj + l — fj + l^O: 


we get, using also (5.8), 

(5.11) (ao(x)v9V(x))/j+i(x) = ^a"(x). 

Thus we have first 


(5.12) 


fj+iix) = 2 


i r 

^ J +00 


<(A) 


(p'+it)aQ{t) 


dt = 
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and, differentiating (5.11), we obtain by induction 


(5.13) 


= 0{{x) 




)■ 


Then (5.5) follows by differentiating (5.10) and using Leibniz formula. □ 


To the formal series Ylj>o defined in Proposition 
a function a by means of a Borel construction, setting 

(5.14) a{x,E,£,h) = aj{x, E, g)/r^x(Aj7i) 

j>o 

for some plateau function x £ over {0}, and a suitable sequence (Aj) 

of real numbers such that Xj —)■ +oo as j —)• +00 (see e.g. m Chapter 2]). 
Then, for all G and any k > 0, 


5.1 


we can associate 


N-l 


(5.15) 




{x,h) — ^ a^^\x)h^\ = 0{h^{x) 




j=0 


Moreover u^^^ikb = is an approximate solution to the Schrodinger 

equation (4.1), in the sense that 

(5.16) Ph,eU+,^kbix, E, e, h) = E, e, h), 


where, for all A^ G N, 

(5.17) r{x,E,£,h) = 0{h^ {x)~^). 


Now we build a solution that has the formal WKB solution constructed 
above as an asymptotic expansion in ]ar + (5o,+oo[. To do so, we estab¬ 
lish first some estimates for the solutions of the inhomogeneous Schrodinger 
equation 


(5.18) 


—h'^u” + (14 — E)u = u. 


on intervals of the form I\ 
write 


A 3A 
^ 2 ’ 


, for large A. For simpler notations we 


(5.19) 


Q{x) = QE,e{x) = Ve-E, 
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and we set 


(5.20) 


X = A + Ax, 


so that X G Multiplying also by A the equation (5.18) on I\ is 

equivalent to the equation 


(5.21) (h^Dl + Q{x))u = V, 
on [—^, 5 ], where 

(5.22) h = Q{x) = A“'”°Q(x), v{x) = A“”^°n(x). 


Notice in particular that Q 1, and that = 0(1) for all k > 1. As in 


Chapter 7], we write (5.21) as the first order system, 
(5.23) 


{hDi + A{x))U = V, A(x) = 


0 

Q(x) 


-1 

0 


[/ = 


u 

hDfU 


V = 


and we denote £ = (£(ij)) G C°°(I x I,A42(C)) the fundamental solution of 
this system. One can prove that (see [23l Theorem 7.1.3]), for any j, k G N, 


(5.24) \\{hdxy{hdx)^£{x,y)\\ < C'j,fcexp(l| Im(^(x) - Im(,5(y)|), 

h 


for some Cj^k > 0 , where is the solution of the eikonal equation associated 
to (5.21) such that is decaying as x increases. Now the solution u 

of (5.21) satisfies, for any x^y G 


(5.25) 


{ t((x) A 

\hDs;u{x)J 


£{x,y) 


] 

\hDiu{y)J 


Since i(p{x)/h 
associated to 
u{x) = 'u(x). 


= i(p{x)/h, where (p is the solution of the eikonal equation 
(5.18|) which decays when x increases, we also have, with 


(5.26) 


u{x) 


= £{x,y) 


u{y) 

jrriQ DyU{y) 


with 

(5.27) \\£{x,y)\\ < C'exp(^| Im(^(x) - lmip{y)\. 

h 
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Replacing A by |x| we get 

where 


(5.29) 




and E satisfies the estimate (5.27) possibly with another constant (7 > 0. 


For general ar + do < x < y,we can cover [x, y] with 1 + C>(ln(|)) intervals 
of the type Ix. Thus, writing 


£{x,y) = £{x,xi)£{xi,x 2 ) ■ ..£{xn,y), 


where each of the intervals [xjjXj+i] is contained in one of the chosen 
we obtain 

(5.30) ||T(x,y)|| < exp(i| Im(^(x) — Im</5(7/)|?/“™'°/^. 

h 

In particular we have, for all ar + 5o < x <y, 

0 ( 1 ) uC’-™o/2 1 

(5.31) |Ti 2 (a;, y)\ = - ~ 

where we have used the fact that for x <y, lui(p(y) > lm(p{x). 


Now we consider again the function u+^wkb, and we denote ux the unique 
solution of the Schrodinger equation (4.1) such that 


(5.32) 


ua(A, E, e, h) = u+,.u,fcb(A, E, e, h), 
hDux(X, E, e, h) = hDu^^^ikbi'^, ^)' 


We have 


(5.33) '^/i,e(l]or+^o A] ("^A '*^+,uifcfe)) l]ar+<5o,A]^) 

so that, for x G]ar + <5o, A], 


(5.34) ux{x) - u+^y,kb(x) = -- 


£(i, 2 ){x,y)r{y)dy. 
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Then, noticing that we can take C > 0 arbitrarily large in (5.31) , we get, 
for all G N and all /c G N, 


(5.35) uf (rc) - 

and, for Ai < A2 G M"*" large enough, 

(5.36) - uf^{x) = 

Thus, the family {u\) converges to some function u+, which is an exact 


solution to ( |4.l[ ), and ( |5.35[ ) gives that, for all G N and all A: G N, 

(5.37) uf{x) - u%Jx) = 0(A"{i)-«e-‘“^(“^)'''-). 


We have proved the main part of the 


Proposition 5.2 Let So > 0, and £0 > 0 be small enough. Let u = 
Qi-v+lh Qjhl be a formal WKB solution in jor + (5o) +oo[ satisfying ( 5.1 ) 
and ( |5.4p . Then, for any {E,£) G D{Eo,£o) x iA(0,eo), the equation (Tfl) 
has a unique solution on]ar + Sq, +00 [ such that, 

(5.38) u+{x, E, £, h) = a+(x, E, e, 


with 

(5.39) 


a+(x, E, £, h) E {x, E, £)h^, 

j=0 


in the sense that 

\/N G N*,VA: G N, 3 CN,k > 0 such that 

N-l 

(5.40) |5^(a+(x,A) - ^ a+(x)/i^)| < . 

j=0 

Moreover u+ belongs to L'^Qar + Sq, +00[, and it is analytic with respect to 
(E,£} g D(Eo,£o} X iA( 0 ,eo)- 


Proof: It only remains to prove that the solution belongs to L^Qor + 
(Io,+oo[). Assumption (A2) implies that there exists C > 0 such that 
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Re(T4(a:) — E) > -^ for all x large enough. Thus we have, for x ^]a 2 , +oo[ 
large enough, 


(5.41) Re(i(/9+(x)) = — Re / {Vs{t) — E)^ dt <——+ C. 


On the other hand, the estimate (5.40) for iV = 0, A: = 0, gives a+(x, E, e) = 
so that M+ G L‘^{]ar + (5o) +oo[). □ 


It is clear that we have the same result for the existence of a solution 
U- G that has, for any do > 0, a WKB asymptotic expansion in 

] — oo, ai — do[ of the form 

(5.42) u-{x,E,e,h) = a-{x, 


where the phase ip- is defined by 

(5.43) i(p-{x,E,£) = f {Ve{t) — E)jdt 

Jai{E,e) 

1 

where we recall that the determination of 11—(14 (t) — E)^ is fixed in such 
a way that, for x g] — oo, — do[. 


(5.44) 


Iie{i(p-{x,E,£)) < 0. 


We also have, in the same sense as in (5.40), 


(5.45) a_(x, E, e, h) ~ a - (x, E, e)h^, 

j=0 


where u^kb = ^ formal WKB solution. It is of course 

also analytic with respect to {E,e) G D{Eo,£o) x T)(0,eo)- 


6 The quantization condition 


To start with, we derive the quantization condition, using only the double 
well structure but not yet the PT-symmetry nor the symmetry following 
from the fact that Pq is real and self-adjoint. 
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In Section]^ we have introduced the two null solutions ri+, U- of [P^ — 
E)u = 0 that decay exponentially near +oo and —oo respectively, and we 
know that E is an eigenvalue of Pg precisely when W (u+, m_ ) = 0 or equiv¬ 
alently when and U- are colinear. It is clear that we can choose u± and 
Uq, Uq of the preceding section, so that 

Uq = U-, Uq = 


( 6 . 1 ) 


By (4.13), (4.22) we have 
4 = 


Similarly, by ( |4.23 [), (|4.24|), 

( 6 . 2 ) 


r 

Uq 




Here, we recall (4.26), (4.28), implying 
(6.3) 


r = (lA — 


Vn = e 


and (4.29), (4.30), that give 
(6.4) 


2i 


*’» = (-1 - ”-i) 


ui - 


(4.31) and (6.3) form a system that allows to express in terms of Uq, Ug. 
Similarly, (4.32) and (6.4) allow us to express u!(_ in terms of Uq, Uq. After 


some straightforward calculations, we get. 


(6.5) 


1 - 


1 (l \ f Vq 


v'_ij 7(1+7^ Vl 


li) Wn ' ’ 


( 6 . 6 ) 


u[ \ _ 1 / 2ie-^l'^Y+ 1 

A J ~ Y-+1+ 1 


Combining (6.1) and (6.5), we get after a straightforward calculation, 
(6.7) 


Uq = 


1 


/ I / 

7 - + 7 + 
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Combining (6.2), (6.6), we get 

( 6 . 8 ) 

1 


r 

Uq = 


l- + l\ 




7-) 



Since Vq, Vq are linearly independent, we see that E is an eigenvalue of 
Pe precisely when the two row matrices in (6.7) and (6.8) are colinear or 
equivalently when the determinant of the matrix, formed by these two rows, 
is equal to 0. We then get the quantization condition 

0 =1 1 {^Irlh _ ^-ilrlh'^ 

which we rewrite as 


(6.9) 

where 


( 6 . 10 ) 


f{E,e) = 0, 


f{E,e) 


— sin (li/h) sin (Ir/h). 


We shall now take into account the various symmetry properties. For 
functions u(x,E,e), where x,E vary in some domains in C and e in some 
real domain, we put, 

(6.11) Co {u){x, E, e) = u*{x, E, e) = u{x, E, —e), 


(6.12) Pt (n)(x, E, e) = u^{x, E, e) = u{—x, E, e), 

so that Pt is equal to VT in the introduction. Notice that Pt and Co are 
idempotent anti-linear operators that commute: Co o Pt = Pt o Co. 

Using only that 14 = Vq+ieW with Vq, W real-valued on the real domain, 
we see that 

v;ix) = = 14(x), Pt (K)(x) = = Vo{-x) - ieWi-x), 
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and if we make the PT-symmetry assumption (A5), we get Pt(f4) = 14 and 
hence that 


(6.13) {Pe-E)oPt =Pt o{P^-E). 

Without (A5), we still have 


(6.14) {Pe-E)oCo = Coo {P^-E). 


To verify this quickly, we observe that Co (uv) = Co {u)Co (v) for products of 
functions and similarly for Pt, and that Pt odx = —dxoPt, Coodx = c^^oCo . 
Recall that to leading order, 


4 = := {Ve-E),-^e-^J-^ 




By straightforward calculations, oberving that a£{E,—e) = a£{E,e) (i.e. 
= ag), we obtain. 


(6.15) 


Co (4,o) = 


i 

0 , 0 - 


In view of (6.14), we know that Co(uq) is a nnll solution oi Pf. — E and 
using also (6.15), we conclude that ^ (wq + Co (uq)) is a null solution with 


leading asymptotics Uq q which is invariant under Co, so if we replace Uq by 
this function we gain the property. 


(6.16) 


Co (4) = 


Similarly, in the discussion leading to (4.12) we see that we can choose 


so that 

(6.17) 


uf = Co (till). 


Since form a basis for the space of null solutions of — 77, we get 

from (4.12)^R15), that Co(r+) = t_, so afte r rep lacing by we 

still have (6.17) for the new functions in (4.13). 


Next, notice that Ig in (4.19) satisfies 

/; = ig. 


vf = Co (nil), 


(6.18) 

This means that 

(6.19) 
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for in (4.22) 


The whole discussion so far applies with ‘T” replaced by “r” and we get 

(6.20) Co (mS) = uS, i; = Ir, u[ = Co (uLi), vl = Co (uCi). 
Moreover, we check that 

(6.21) J* = J. 


From (4.26), (4.28) we now get 

( 6 . 22 ) 


Similarly, 

(6.23) 


Co(uS) = < 


Co (u^) = vl 


Then in (4.31), (4.32) we must have 
(6.24) Co( 7 ^)= 7 i, Co ( 7 ;) 


= 7-- 


It follows that / in (6.9), (6.10) satisfies 
( 6 . 25 ) r = /• 


Let us finally use the "PT-symmetry assumption (A5) or equivalently 


(6.13). We then check that 


(6.26) 4 = Ir, 
and 

(6.27) = J, 
and that we can choose 

(6.28) Uq = Pt (ug), = Pt (uli), = Pt (uf), 


(6.29) 


4i = 
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Then by (6.3), (6.4), 
(6.30) 


ng = Pt (n^). 


and from (4.31), (4.32) we infer that 

7i=Pt(7i). 


(6.31) 

It follows that 

(6.32) 


P = f- 


We have seen that the zeros of coincide with the eigenvalues of Ps in 

a neighborhood of Eq . We end this section by showing that the multiplicities 
agree also. 


Recall (6.7), (6.8) that we write as 


Uq = a^VQ + b^v^, 


Taking the Wronskians, we get 

Here the last Wronskian is non-vanishing, so up to a non-vanishing holomor- 
phic factor / is equal to W{uq,Uq) and the zeros of /, counted with their 
multiplicity coincide with those of W{vq,Vq). 


Hence it remains to identify eigenvalues of counted with their multi¬ 
plicity with the zeros of the Wronskian W{uq,Uq). For that we can widen 
the perspective slightly and apply a general discussion: 


Let o G M. Let A vary in neigh {Eq, C). (The symbol “17” will temporarily 
be used to denote operators.) Using the ellipticity of Pg — A near -l-oo, we 
see that the right Dirichlet problem 

{Pe — X)u = V, u{a) = u+ 

has a unique solution u G H‘^{]a, -|-oo[) for every {v, u_|_) G H^{]a, -|-oo[) x C. 
Similarly, by using the ellipticity near —oo we see that the corresponding 
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left Dirichlet problem has a unique solution u G H‘^{] — oo,a[) for every 
— oo,a[) X C. Denote the solutions to the two problems by 
u = ErV + v+ and u = Egv + E^ respectively. 

It follows that the Grushin problem 

f {Ps — X)u + R-U- = V, 

= v+, 


has a unique solution 

{u, U-) G — oo, a[) © H^{]a, +oo[)) n x C 

for every (u,u_|_) G H^(R) x C, where 

R+u := u{a), R-U- = uSa-, 

and 5a denotes the delta function at x = a. Indeed, the solution is given by 

{EiV + Ejv+, X <a 2- , n^ 

u{x) = < u- = [h Ozu){a — 0) — [h dzu)[a + 0 ). 

yErV + Ejv+, X > a, 

We write this solution, 

u\^(E E+\(v\ 

u.) \E. E.+J \v+) • 

It is a standard fact for the Grushin reduction (see m Section 6, Appendix 
A], or [23]) that the eigenvalues counted with their multiplicity, coincide 

(near Eq) with the zeros of E _ ^(z) counted with their multiplicity (e is here 

fixed and suppressed from the notation most of the time). For completeness, 
we recall the proof. Let zq be an eigenvalue of P. Its multiplicity m{zo) is 
equal to the rank, and hence to the trace of the spectral projection 

m{zo) = tr ^ / (^ - Py^dz, 

where 7 is the oriented boundary of a small disc centered at zq. Now 
(z - P)-^ = -Eiz) + E+{z)EZliz)Eyz), 
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and E{z) is holomorphic near zq, so 


(6.33) 


m{zo)=iT^ j E+{z)Ezl{z)E_{z)dz 
= ^ / ty:iE+{z)Ezliz)E_{z))dz 
= ^ / ty:iEzl{z)E_{z)E+{z))dz 
= ^ / Ezl{z)E_{z)E+{z)dz. 


Finally, since E-{z)E^{z) = dzE_^{z), we see that m{zo) is equal to the 
multiplicity of zq as a zero of E _ 

Denoting u_ = £'^(1), u+ = £+(1), we notice that 

£_^ = h?{dzU-{a) — dzUj^{a)) = /iPF(«_,u+), 

since u^{a) = 1. Thus, in a neighborhood of £o, the eigenvalues of Pe 
counted with their multiplicity can be identihed with the zeros of A i—)■ 
lT(n_(A),u+(A)). 


If neigh (£o, C) 9 A i—)• Uzp{z, A) are holomorphic families of null solutions 
to ££ —A, exponentially decaying near =Foo and ^ 0, VA, then Szp = cJzp(A)uqi, 
where Uip are holomorphic in A and non-vanishing. Thus the zeros of 
W (rt_, tt+) = cj_(A)cj+(A)IT(u_, u+) coincides with those of W (n_, u+) and 
we have completed the identification. 


7 The behaviour of the eigenvalues 


In order to study the zeros of/(•, e), we hrst recollect the various symmetries: 


(7.1) 

II 

II 

* 

II 

j* 

±r 

— I — I 

and 





(7.2) 

1 

* 

II 

-J 

T)* = 

V 

7+> 

II 
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where we recall that 7 ^ = 1 + 0{h). 

Let us first look at the factor 

(7.3) g = ge{E, s) = 


and drop the super/subscript i when convenient to do so. From (7.1), (7.2), 
we infer that 


(7.4) 
Write 

(7.5) 


9 =9 


1 1 
..i \ 2 / i \ 2 


7i=(7i7;j (TiMj 
7^=^7l7^)'f7l/7^)■^=:p.e-^^ 


where we choose the branches of the square roots close to 1 and the logarithm 
close to 0, so that p = 1 + 0{h), 9 = 0{h). Then 

(7.6) p* = p, 
so (e*®)* = and hence 

(7.7) 9* = 9. 

We write 

(7.8) p^ = p^(eK^^ + e-s^0, h = h + h9e = h + 0{h^) = 1}. 


Similarly, we consider 

(7.9) g = gr(E, e) = (eS^’'7+ + • 


and again we have (7.4), now with g = g^. Write 
(7.10) 


7; = (7!:7;)M7;/7d)^=: 


again with p = 1 + 0{h), 9 = 0{h) satisfying (7.6), (7.7). 
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We write 


(7.11) g^ = Ir = Ir + h 6 r = Ir + 0 {h^) = i;. 

We now take into account the VT symmetry. Clearly, 

(7.12) pI = Pr 

and from (7!)^ = 7-) we get pre~^^i = so 

(7.13) 9 ^ = Or and hence = Ir- 


Using also that jj = Ir, we can rewrite / in 

/(E.e) (ei' + e-if) (eU' +e-i'') 

- ie-2'^/'^sin(///i)sin(/t/h) 

=pp^ cos{I/h) cos{I^/h) — -e sm{I/h)sm{P/h) 

where I = I = Ii, p = p£. Dividing this function with ppl will not 
modify the zeros and we get the new (slightly modified) function that we 
shall denote by the same symbol, 

(7.14) f{E,e) = cos(///i) cos(/Vh) — sin(///i) sin(/'l’//i), 

where 

(7.15) J = J + /iln(pp+) = J + C>(/i2). 

We know that 


(6.10) as 


(7.16) 
From 

(7.17) 
we get, 

(7.18) 


P = f = r- 

r0i 

I{E,e)= / {E-Vs{x))'2dx, 

Jai 

dEl{E,e) = l^f {E-Ve{x))~^dx, 

^ J Oil 
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and 


1 


deIiE,e) = - / iE-V,ix))—2W{x)dx. 
2 * Ja, 


(7.19) 

It follows that 

(7.20) dEl{E,0) > 0, icle/(i7,0) G M when E G neigh (£'o,] 


We now adopt the assumption (A7) so that the integral in (7.19) is non¬ 
vanishing for {E,s) = {Eq,0) and in order to fix the ideas (possibly after 
replacing (e, bh) by (—e, —W)) that 


(7.21) 

so that 

(7.22) 


/ (So-yo(x))-5W(x)dx>0, 

J OL^ 


idsI{E, 0) > 0 for 77 G neigh {Eq,R). 


Let Jo = I{Eq, 0) G M. In view of the first part of ( 7.20[ ), the map ■ 

neigh {Eq, C) —)• neigh (Jq, C) is bijective for e G neigh (0, M) with an inverse 
K{-,e) : neigh (Io,C) —)• neigh (i7o,C) such that K{i,£) is holomorphic in 
(i, e). We also know that K{i, 0) is real when l is real. The property I* = I 
implies that I{E,e) is real when e = 0 and since 1 = 1 + 0{h'^), we see 
that 7(-,e) has a local inverse K with the same properties as K. Further, 
K = K + 0{h?) has a complete asymptotic expansion in powers of h in the 
space of holomorphic functions defined in a neighborhood of (lojO). 


The zeros of the factor cos{I/h) are given by 


(7.23) 


I{E,e) 



nh, 


for /c G Z such that [k + irh belongs to a neighborhood of Iq. (The 
classical action for the left potential well is equal to 21, so we recognize the 
Bohr-Sommerfeld quantization condition 2/(77, e) = {k + l/2)27r/i.) They 
are situated on the real-analytic curve 


(7.24) r(e) = {Ee neigh [Eq, C); I{E, e) G M}. 


Alternatively, the zeros are of the form 

(7.25) Ek = K{{k + l/2)TTh,£), 
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and the curve (7.24) is of the form r(e) = iir(neigh (0, M), e). 


The curve T can also be represented in the form 

T: ImE = g(ReE,e), 

where 

9{t, e) ~ git, e) + hgiit, s) + 

and 

r : Im£' = g(ReE,s) 

is the curve, determined by the condition I{E,£) G M. We know that T, T 
are real segments when e = 0, so 

g, g = 0 {£). 


Writing 

Im/(Rei7 + ig(Re E, £), £) = 0, 
and differentiating with respect to e at e = 0, we get 

7 f) I 

dsgiReE,0) = ^(ReE,0)>0. 
OeI 


Hence, by Taylor expansion, 


idj. 


g{ReE, e) = o)e + C>(e^). 

OeI 


It follows that 

g{ReE,£) = 


idel 

OeI 


{ReE,0) + Oih^) e + 0(e2). 


Similarly, by differentiating the equation /(i7fc(e),e) = (k + lj2)'Kh with 
respect to £, we get 

deEki£) = , 

OeI ^ ^ 

and again by Taylor expansion at e = 0, 


Eki^) 

(7.26) 


Eki0) + i'M(Eki0),0)£ + Oi£^) 

OeI ^ ^ 

(EkiO) + 0{h^)) + i (^fc(O), O) + Oih‘^)^ £ + 0(£^) 
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Here Ek{0) and -Efc(O) are real and given by the Bohr-Sommerfeld conditions 
7(^fc(0), 0) = {k + l/2)7r/i, /(^fc(0), 0) = {k + l/2)7r/i. 


The zeros of the second factor 

cos{P/h){E,e) = cos{I/h){E,e), 

are given by -E = Ek{e) and they are situated on the complex conjugate 
curve 

ImE = —^(ReE, e). 

We next make an exponential localization of the zeros of /. We have, 

I coszl X min(dist {z, (Z + l/2)7r), 
so for E G neigh (Eq, C), 

|cos(7/h)| xminQdist(E,{Efc}),l^ 

Since 7=/ + 0(/i2), we get 

(7.27) |cos(7/h)| xminQdist(E,{Efc}),l^ 

Similarly, 

(7.28) |cos(7t//i)| xminQdist(E,{Efc}),l^ 

(7.29) |sin//h| < 2el^“^l/7 
and 

(7.30) |sin/t/h| < 

We conclude that /(E, e) / 0 when 

I cos(7/h)cos(7t//i)| > sin(//h)sin(lt/h)| 
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and that this holds when 


(7.31) _ 
min Qdist (£', {Ek}), min Qdist (£', {Ek}), 

Thus if C > 0 is large enough, then f{E,e) ^ 0 when 

(7.32) E^[jD{Ek, Che- EE,e)/h^ ^ J KeJ{E,e)/hy 

k k 


Now we observe that for any E G neigh {Eq, C), 


After doubling the constant in (7.32), we conclude that 

(7.33) _ 


If .B G T» (^Ek, and if 


(7.34) 


e >/ie 


then dist {E, {En}) e, and from (7.31) we conclude that f{E, e) / 0 if 
(7.35) idist(E,.Bfc)min(|,l) > e-2Rej(Tfe,£)/h^ 


Thus, the zeros of / in D{Ek,Che Rej(£;fe)/h^ contained i 


m 


E Ek, — 


Ch 


^-2Re J{Ei.,e)/h 


min(e//i, 1) 


If we drop the assumption (7.34), we have 


(7.36) 


f-\0,s)f^D (^Ek,Che-^^'^^^'^E/h 

C D {Ek,Chmhi fl,max(V£,l)e"^'^^^-®"’^^/^) ^-KeJ(E^,e)/h \_ 
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The same discussion is valid with Ek replaced by Ek and we get the following 
improvement of (7.33|): 


(7.37) r\^.e) c[jD [Eu,r{Ek,e)) (Eu,r {Ek,e)) , 


where 


(7.38) r{E,e) = Chmin (l,max(Ve, 


By inserting a deformation parameter 6 G [0,1] in front of the second 


term in the last expression for / in (7.14), we will not change the local 


ization (7.37) of the zeros and the number of such zeros in each connected 


component of the set in the right hand side of that inclusion is independent 
of 0. It follows that 


— when these discs are disjoint, /(•,£) has precisely one zero in each of 
D (Ek,r and D (Ek,r (Ek,e] 


— in general /(•,£) has precisely 2 zeros in 

D (Ek,r (Ek,e)] U D (Ek,r [Ek,£ 


From the first of these observations, (|7.33|), (|7.26|) and the fact that 

(7.39) 
we get 


id I 

—^(Fi, 0) > 0 when E is real, 
OeI 


Proposition 7.1 Assume (7.21) for E G neigh (Flo, K) and that e is real 
and 

1 > |e| > 

for some positive constant C. Then the eigenvalues in neigh (Flo, C) are 
simple and non-real of the form Zk{£] h), Zk{£', h), k G Z, where 

Zk = Ek{£-, h) + 


and we recall (7.26), (7.25). The term 0 (h)e Re j(Ei,,£)//i replaced by 


0 {r{Ek,£)), where r{E,e) is defined in (7.38). 
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It remains to make a more detailed study, when 


e| < /ie-(ReJ(Eo)-i/C)//i^ 


and for that we shall view Ek{0]h) as the nondegenerate local minima of 
fo{E, 0), E G neigh (£'o,M), where we put for 0 < 0 < 1, 


(7.40) 


fe{E,£) = cos{I/h) cos{P/h) 


^e-2^/'*sin(I/h)sin(/t//i), 


so that /i 
real E 


f in (7.14). 


Using that fg 


fe and 4 = fe 


we see that for 


— Jq is real-valued, 

— fo is an even function of e. 

Write 


/o = 9 {E, £)g\E, e), g{E, e) = cos{I/h) = g*{E, e). 

Let Ec{0) be a (real) zero of g{E,0), so that Ec{0) = Ek{0) for some k G Z. 
For e = 0, we have /o = g{E, 0)^ > 0 and Ec{0) is therefore a nondegenerate 
local minimum of /o with /o(i7c(0),0) = 0. Extend Ec{0) to an analytic 
family Ec{e) of critical points of /o(-,e): 

(7.41) dEME,{£),£) = 0. 


We have for real E: 

dEfo{E,s) = 2Re{dEg{E,£)g{E,e)), 
and differentiating this once more and putting e = 0, E = Ec{0), we get 
9|;/o(Ee(0),0) = 2dEg{Ec{0),0)dEg{E,{0),0) = 2 {dEg{Ec{0),0)f , 


i.e. 


(7.42) 


/i'9|/o(^c(0),0) =2 (sin(/(Ee(0),0)/h))' (9i,/(E,(0), 0))' 


=2 


(5eI(E,(0),0)) , 
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where the last identity follows from the fact that cosI{Ec{0),0)/h = 0. 


Differentiating (7.41[), we get 


(5|/o) deEc + dedEfo = 0. 


Here we recall that when E is real, fo{E, e) and 9 e/o(-E, e) are even functions 
of e and hence d^dEfoiE, 0) = 0. It follows that 


(7.43) 


{d,dEfo){E,{0),0) = 0, 


(7.44) d,E,{0) = 0. 

Using this, we get 

ide)l=oifoiEcie),e)) 

= {de),^o f (<9i;/o) {Ecie),e) deEcie) + (defo) {Ec{e),e) 

V =0 

= MEM, 0) + {dEdJo){E,(e),e) deE,{0) 

=0 =0 

= 2\degiE,{0),0)\\ 


Thus, 


(7.45) 


(4)Lo iMEcie),e)) =2 sin/(i?e(0), 


del) (i?c(0),0) 


=2 


dj) (i?c(0),0) 


We next extend Ec{s) to an analytic function Ec{s, 6) determined by the 
conditions Ec{e,0) = Ec{e), 

(7.46) dEfB{Ec{e,e),e) = 0. 


It will be convenient to restrict the attention to a window of size 0{h): 
(7.47) E = Ei + hE, s = he, 
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where Ei G neigh (Eq,M) is a parameter and E G neigh (0, C), e G neigh (0, M 
are rescaled variables. It will also be convenient to have a “Floquet param¬ 
eter” K G M and introduce the following extension of (7.40): 


(7.48) 

fe{E,e,K) = 


cos \ ~ j 


h 


e 


— K I “ sin I ^- K \ sm 


h 


— K 


which coincides with f 0 {E,e), when k G vrZ. Again, fo is real-valued when 
E is real and an even function of e. If we let Ec{£, n, 9) denote a local 
minimum of fg{-,e,K), then (7.42), (7.43), (7.44) extend naturally to the 
case 9 = 0. Writing 

K = K + I {El, 0)/h, 

we get 

(7.49) fg = fe{F,e,K;h) = a{E,£,K-,h) + 9e~‘^EEifl)/hb(^F^e,K-,h), 
where a, b are classical symbols of order 0 in h: 

a ~ ao -|- hai + ..., 5 ~ 6o + hbi + ... 


The critical points with respect to E are nondegenerate and their number 
is uniformly bounded. They have asymptotic expansions in powers of h of 
the form, 

Ec{e,K,9-h) ~ E^{e,K,9e-‘^EEifl)/h)-^FhFl{£,K,9e-‘^EEifi)/h)^ ^ ... 


which gives 

(7.50) Fc{£,K,9;h) = F^{£,K-,h) + 9e-^'^^^^’^'^/’^F^{e,K,9-,h), 


where F^ are classical symbols of order 0 in h and also holomorphic func¬ 
tions. Notice that the terms in the asymptotic expansion of F^ in powers of 
h are independent of 9. Ec = Ei + hFc will be a critical point of fg{-,£) in 
(7.40) when k G vtZ, i.e. when 


(7.51) 


I{Ei,0) 


mod vrZ. 
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It follows from ( |7.50[ ) that 

Fc{e,K,e-,h) = Fc(e,K, 0 ;/i) + 
and hence that 

(7.52) Ecie, K, 9) = Ec{e, k, 0) + 

To evaluate the critical value fe{Ec{e,K,9),£,K) =: fg{£,K) for 0 = 1, we 
notice that 


de{f9i£,K)) = {defe){Eci£,K,9)) 




h 


It 


sin - - K sin — - K {Ec{£, k,6),£, k) 


h 


^ _ :(g-2J(Ec(0A,0),0), 


sin -- K 


{Ec{0,k,0),0,k)) + 0{1) 


£ 

+ V 


Here, we use that 


I 


sin v — (EciO, K, 0), 0, k) = ±1, sin ( v — ^ = sin t ~ ^ ^ 


h 


I 


h 


h 


and integrate from 0 = 0 to 0 = 1, to get 

fii£,K) = /o(e,K) 

(7.53) 1 _ 


_ lg-2J(£;c(0A,0),0)//i _|_ 0 ^g-2Re J((£;i,0),0)//i _|_ M _|_ 0{h)^ ^ . 


Now, return to the window (7.47), where fg 


(7.49) and the critical point Fc is as in (7.50). 
suppressing the corresponding subscript 1) 


= fg{F,£,K;h) is given by 
We have with 9 = 1 (and 


(7.54) f{£,K]h) 


f{£,K;h) 


gi{£,K;h) +g2i£,K-,h)e 


2J{E^{0,K,0)fl)/h 

? 


where gj are classical symbols of order 0 in h. (We first get this with 
J{Ei,0)/h in the exponent, but the replacement by J{Ec{0, k,0),0) does 
not modify the general structure of the formula.) (7.53|) shows that 


(7.55) 


51(0 , k;/ i) = 0, 52(0 , k) 


1 

4’ 
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where gj^ is the leading term in the asymptotic expansion of gj. 


From (7.45), we deduce that 

(7.56) 5|5i,o(0,k) = 2|9,I(F;i,0)|2 > 0. 


Combining ([7^, ([7^, ([7^, we get by Taylor expansion, 

(7.57) J%£, k; h) = 52 ( 0 , k; + k{£, k; h)t^, 


where A: is a symbol of order 0 in h, holomorphic in the other variables, even 
in e and satisfying /c(0,k;0 ) = \d£l{Ei,0)\‘^. has precisely two 

zeros in a neighborhood 0 which are real and of the form E£c{k\ h), where 

(7.58) ec{K-, h) = i{K-, 


and £ is a symbol of order 0 with leading term 

^ 2\d£l{E^,0)\' 

Using that our functions are holomorphic in e, k, we see that 

(7.59) h) = m(e', k; h){^ — £c(ii; h)^), 

where m is holomorphic in e, k and a symbol of order 0 in h with leading 
term niQ (e, k) , satisfying 

(7.60) mo{0,K) = \deI{Ei,0)\. 


( |7.42 ) can be extended to the K-dependent case: 

(7.61) {h^dlfo){E£iO, K, 0), K, 0) = 2{dEl{E£{0, 0 ))^. 

This implies that 

(7.62) {dlfi){Ec,e,K;h) = idEl{Ec{0, K,0),0)f + 0{h), 

where the remainder has a complete asymptotic expansion in powers of h. 
By Taylor expansion, 

(7.63) fi{F,e,K;h) = f{£,Ti]h) + g(F, e, k; /i)(F - Fc(e,K, l;/i))^. 


52 















where g > 0 is a symbol of order 0 and 


(7.64) 


q{F,, 0 , 0 ) = 2{dEl{Ec{Q, k, 0),0)f. 


f{-,£,K;h) has two zeros in a small neighborhood of Fc when counted 
with their multiplicity: 

— When |e] < SciE, h) the zeros are real and simple, given by 


(7.65) q{F, £, K] h) 2 (F - Fc(e, k, l]h)) = ±(-/'=(e, k; h)) 2 


When |el = ec(«^; h) we have a double zero, 


(7.66) 


F = F,. 


— When |e] > £c{K',h) the zeros are non-real and simple and complex 
conjugate to each other, given by 
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